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Correlating spin transport and electrode magnetization in a graphene spin 
valve: simultaneous magnetic microscopy and non-local measurements 
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Using simultaneous magnetic force microscopy (MFM) and transport measurements of a graphene spin valve, 
we correlate the non-local spin signal with the magnetization of the device electrodes. The imaged magneti¬ 
zation states corroborate the influence of each electrode within a one-dimensional spin transport model and 
provide evidence linking domain wall pinning to additional features in the transport signal. 


Electrical injection and detection of spin have become 
common techniques to study spin transport and relax¬ 
ation in metals and semiconductors 1 6 and are attrac¬ 
tive for technological implementation of spin-based logic. 
Graphene has demonstrated great promise as the spin 
transport channel in such devices, on account of its rel¬ 
atively long spin diffusion length and lifetime at room 
temperature 4,7,8 . Optimization of these devices necessi¬ 
tates an understanding of the fundamental interactions 
governing spin transport and will require a variety of 
measurement tools and techniques. The successful com¬ 
bination of scanning probe microscopy with charge trans¬ 
port measurements 9-14 suggests great opportunities for 
similar integration with spin transport devices. The 
imaging mode demonstrated in this letter—simultaneous 
MFM and spin transport—provides independent and 
complementary probes of magnetization, and facilitates a 
more complete understanding of spin and magnetization 
coupling in nanoscale devices. This approach provides 
a useful tool for studying complex magnetization con¬ 
figurations and their relationship to transport in other 
technologically relevant devices 15 . 

Here we report direct imaging of electrode magneti¬ 
zation of an operating non-local graphene spin valve. 
These images allow unambiguous determination of the 
overall magnetic state of the device and can be corre¬ 
lated with simultaneously acquired non-local magnetore¬ 
sistance. Using this correlation, we verify that a one¬ 
dimensional spin diffusion model 16 , combined with elec¬ 
trode magnetization switching, quantitatively captures 
the primary behavior of the device—high, low, and in¬ 
termediate resistance states. 

Direct correlation of imaging and transport also pro¬ 
vides insight into the origins of various observed features 
in the transport signal beyond simple binary switching 
of electrode magnetization. We observe domain wall pin¬ 
ning and de-pinning in the ferromagnetic electrodes con- 
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current with otherwise obscure transitions in the non¬ 
local signal. Without imaging, the origins of these 
occasionally-detected transitions 4 would remain unclear. 
Our observation adds to a physical understanding of 
these sources of noise-like variations and aligns with 
the previous observation that ferromagnetic contacts can 
have a strong influence on spins in the transport channel 
via stray field effects 17 . 

The graphene spin valves used in this study consist 
of exfoliated single layer graphene with Co/MgO ferro¬ 
magnetic contacts (see Ref. 18 for details). Using a com¬ 
bination of two-point and four-point measurements, we 
find the particular device described in this manuscript 
exhibits a graphene sheet resistance of 1.38 kD/D and an 
average contact resistance of 6.13 kU. 

Spin transport measurements are performed using the 
traditional four-terminal non-local scheme 3,4 (see device 
schematic in Fig. 1(a)). All measurements are performed 
at room temperature and in vacuum. We source 10 pA 
of current at J + (drain at /_) at 11 Hz, and measure the 
non-local voltage Unl = V+ — V- with a lock-in amplifier. 
Results are plotted as Rnl = Unl/U An external mag¬ 
netic field is swept in-plane, parallel to the long axis of 
the Co electrodes, in order to obtain the non-local mag¬ 
netoresistance (NLMR) signal (Fig. 1(b)). From a Hanle 
measurement of this device (with magnetic field oriented 
perpendicular to both the substrate plane and the in¬ 
jected spin orientation), we obtained a spin lifetime r s = 
476 ps and spin diffusion length A s = 3.10 pm by using 
the measured device resistances and the fit algorithm of 
Ref. 19. These parameters are used as inputs to the spin 
diffusion model to produce the spin electrochemical po¬ 
tential (ECP) profiles shown in Figs. 2, 3, and 4 and to 
fit our NLMR data (Fig. 1, solid lines). 

We have previously described the custom-built scan¬ 
ning probe microscope which enables simultaneous acqui¬ 
sition of force microscopy and transport measurements 20 . 
This system utilizes frequency-shift detection of an os¬ 
cillating cantilever to image atomic, electrostatic, and 
magnetic forces. The magnetic force microscopy (MFM) 
signature of a bar magnet (which adequately approxi- 
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mates the Co electrodes) creates opposite frequency shifts 
at opposing ends of the bar, where the field gradient is 
strongest. This MFM fingerprint enables clear determi¬ 
nation of the magnetization orientation of each electrode. 
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FIG. 1. (a) Experimental schematic showing non-local cir¬ 
cuit configuration of a graphene spin valve and the scanned 
magnetic probe positioned above the device. The device con¬ 
tains 7 magnetic Co electrodes, 4 of which are used for trans¬ 
port measurements. One possible magnetic configuration is 
shown with colored arrows representing magnetization direc¬ 
tion. (b) Electrically-detected non-local spin valve signal. At 
magnetic field values indicated by numbered triangles, we ac¬ 
quired MFM images. A fit to the data is performed using the 
MFM-determined magnetic configuration of each state (solid 
lines). 


At specific fields during the NLMR measurement field 
sweep (indicated by numbered triangles in Fig. 1(b)), 
we acquired MFM images of the graphene flake and the 
nearby ends of all seven Co electrodes. Figs. 2 and 3 
show the sequence of MFM images, constituting direct 
observation of magnetization orientation and switching 
sequence in an operational non-local spin valve. These 
magnetization states can be explicitly correlated with the 
measured NLMR. 

To begin the measurement sequence, the field is first 
ramped to —250 mT to initialize all electrodes to the par¬ 
allel state. The up sweep non-local voltage data is then 
recorded starting at — 50 mT. The first MFM image is 
acquired on the up sweep at 0 mT applied field, verifying 
the all-parallel configuration (Fig. 2, subpanel 1). The 
graphene flake is outlined in dashed lines for reference. 
It is faintly observable in the force detection images ow¬ 


ing to the periodic 11 Hz voltage applied to drive the 
injection current. This sinusoidal voltage is experienced 
by the scanned cantilever as a periodic frequency shift 
due to the capacitive interaction between the grounded, 
conducting cantilever and the biased sample 21 . 
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FIG. 2. LEFT: MFM images acquired during magnetic field 
up sweep. Color scale represents measured cantilever fre¬ 
quency shift, caused by tip-sample forces. Images corre¬ 
late with the Rnl measurements in Fig. 1(b) indicated by 
black numbered triangles. RIGHT: 1-D spin diffusion cal¬ 
culation, showing spin-resolved electrochemical potential as 
a function of position in the graphene channel for up (blue) 
and down (red) spins. Positions of the four circuit electrodes 
(V_,V + ,J+, and VI) are indicated with dashed lines, and the 
overall device state (magnetization direction of electrodes) is 
indicated with up (blue) and down (red) arrows in the in¬ 
set. The measured voltage difference Vnl is proportional to 
the spin density difference V+ — V- highlighted in these 1-D 
calculations. 

On the up sweep, 5 distinct magnetic configurations 
are observed via MFM, corresponding to four unique Rnl 
levels (the all-parallel configurations in subpanels 1 and 
5 exhibit the same resistance). From this sequence of 
images, and correlation with the measured Rnl levels 
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at the corresponding field, it is clear that all four circuit 
electrodes play critical roles in determining the measured 
non-local voltage. The injector (/ + ) and extractor (/_) 
contacts determine the steady-state spin polarization in 
the channel, while the detector (U + ) and reference (V_) 
contacts probe the spin chemical potential set by their 
respective magnetic orientations. 
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FIG. 3. Sequence of MFM images and spin ECP profiles 
correlated with Rnl measurements at specific field values (red 
numbered triangles in Fig. 1(b)) acquired on magnetic field 
down sweep. The MFM tip magnetization has reversed prior 
to imaging state 6, inverting the frequency shift color contrast 
relative to Fig. 2. The electrode magnetizations are indicated 
by colored arrows in the spin density profile insets using the 
convention established in Fig. 2. 

A similar set of states are observed on the down sweep, 
as shown in Fig. 3. The electrode magnetizations are 
again initialized to an all-parallel state with a large 
+250 mT field. This field reverses the MFM tip mag¬ 
netization, which causes an inversion of the frequency 
shift contrast in Fig. 3 relative to Fig. 2. Therefore, the 
electrode magnetizations in states 5 and 6 are identical. 

With knowledge of the total magnetization state de¬ 
termined by the MFM images, we can model the spin 
ECP in the channel to understand the magnitude and 
sign of the measured non-local voltage. The spin ECP 
spatial profile for each of the 9 observed electrode config¬ 
urations is calculated using a 1-D spin diffusion model 16 . 
These are shown immediately adjacent to the MFM im¬ 


ages in Figs. 2 and 3 in order to highlight the correla¬ 
tion between the imaged magnetization orientation and 
the corresponding spin profile in the graphene. Inputs 
to the model include the geometry of our device, the 
switching fields of each electrode as determined by MFM 
and transport, and the spin diffusion length (determined 
from a Hanle measurement). The 1-D model does not 
include contact-induced spin relaxation 22 (i.e. once in¬ 
jected, graphene spins cannot escape through the metal¬ 
lic electrodes). We find this to be justifiable given mea¬ 
sured device resistances (^contact > /^graphene) an d the 
good agreement between experiment and fit. Using this 
model, we can replicate the 6 unique voltage levels ob¬ 
tained in the electrically-detected non-local spin signal as 
shown from the fit in Fig. 1 (solid lines). There are no 
free parameters in the fit aside from a scaling of the over¬ 
all magnitude, and a DC offset to match the measured 
data. The relative voltage levels of the various resistance 
states are fixed by the experimentally-determined model 
inputs discussed above (i.e. the voltage level of each state 
cannot float freely with respect to the others). 

The MFM images also uncover device behavior beyond 
the model of single domain switching. We can correlate 
several transitions in the NLMR signal (Fig. 4(a) and 
(d), green triangles) with interesting features in MFM 
images obtained at the same field. This direct imag¬ 
ing provides evidence linking magnetization domain wall 
pinning and de-pinning to the observed changes in the 
NLMR signal. Uncovering this connection between do¬ 
main wall motion and spin transport is a key new capabil¬ 
ity made possible by simultaneous imaging and correlated 
transport. 

Fig. 4(b) shows a domain wall pinned directly above 
the graphene channel on the Co electrode between U_ 
and U + at a field of — 13.6mT. This domain wall pro¬ 
duces a large, but localized, stray field experienced by 
spins in the underlying graphene. These spins are de- 
phased by the field, causing a change in Rnl indicated 
by the green triangle in Fig. 4(a). Following the meth¬ 
ods of Ref. 23, we can model spin transport in the pres¬ 
ence of such a localized, inhomogeneous field to obtain 
the new steady-state spin density. Qualitatively, we re¬ 
produce the observed reduction in spin signal. We find 
that the dephasing has a greater impact on spins be¬ 
low U + than those under U_, thereby reducing the volt¬ 
age difference Vnl (see spin density profiles, Fig. 4(c)). 
Growth of Co on graphene can result in enhanced grain 
and domain formation of Co due to graphene’s high sur¬ 
face diffusion 18,24,25 , and domains could be pinned above 
the graphene by these grain boundaries. 

A similar feature is observed during the up sweep in 
Fig. 4(d), between 8 and 9.4 mT. At 9.4 mT, we acquired 
the MFM image in (e) showing that the leftmost elec¬ 
trode (external to the non-local circuit) had fully reversed 
relative to the OmT image of Fig. 2, subpanel 1. Prior 
to this reversal, a domain wall necessarily propagated 
across the graphene channel which would have produced 
a stray field experienced by the graphene spins. Modeling 
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FIG. 4. (a) Down-sweep NLMR signal. Green triangle in¬ 
dicates field value at which the MFM image shown in (b) 
was acquired. In this MFM image, a pinned domain wall 
can be clearly seen on the electrode between V- and V+. (c) 
Modeled spin ECP with (Fnl,dw, green) and without (Vnl, 
black ) the pinned domain wall. The stray field from the 
domain wall causes local dephasing and a reduction in Vnl 
(green curve), (d) Up-sweep NLMR signal. Green triangle 
indicates field value at which the MFM image shown in (e) 
was acquired. In this MFM image, the outermost electrode 
has reversed (relative to the state 1). Immediately before this 
electrode’s magnetization reversed, the NLMR signal was seen 
to increase, and then return to baseline upon reversal. 


spin diffusion in the presence of this localized field repro¬ 
duces the qualitative behavior of the electrically-detected 
NLMR: the non-local voltage is enhanced by the dephas¬ 
ing effect of the stray field (since spin accumulation near 
V- is primarily affected). This enhancement persists un¬ 
til the electrode has completely reversed at 9.4 mT and 
the stray field has been eliminated. 

In conclusion, we have utilized simultaneous transport 
and magnetic force microscopy to correlate the com¬ 
plete magnetization state of a non-local graphene spin 
valve with its spin signal. We find good agreement be¬ 
tween measured voltages and a 1-D spin diffusion model 
informed by the MFM images. Such direct correla¬ 
tion studies also uncover device behavior that would be 
impossible to understand with transport measurements 
alone—namely, spin transport sensitivity to pinning and 
de-pinning of ferromagnetic domains in the contact elec¬ 
trodes. 

These results point to possible studies and applica¬ 
tions in which spins interact with a much richer set of 
states than the two allowed orientations of a monolithi- 
cally magnetized ferromagnet. Understanding the inter¬ 
actions between spin transport and mobile domain walls 
or other magnetic textures in nanoscale geometries is an 
area of intense current research 15 . In-operando MFM of 


functioning devices is well suited for such studies. For 
example, the fast and sensitive force detection provided 
by our microscope 20 is able to detect domain wall tran¬ 
sit below the scanned tip, which could be correlated with 
high-bandwidth features in transport data. Furthermore, 
the local spin dephasing due to stray fields from the con¬ 
tact electrodes is identical to the spin imaging mecha¬ 
nism used in scanning spin precession microscopy 23,26 , 
and suggests that spins in graphene should be amenable 
to such an imaging technique. Overall, this approach 
shows clear promise for the precise characterization of 
spintronic device performance, and the development of 
reliable, application-ready devices. 
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